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We study the performance of achromatic half-wave plates (AHWP) as a 
function of their construction parameters, the detection bandwidth of a power 
detector operating in the millimeter wave band, and the spectral shape of 
the incident radiation. We focus particular attention on the extraction of the 
degree of incident polarization and its orientation angle from the intensity 
measured as a function of AHWP rotation angle, which we call the IVA (in- 
tensity versus angle). We quantify the phase offset of the IVA and point to 
potential systematic errors in the extraction of this offset in cases where the 
incident spectrum is not sufficiently well known. We show how the phase offset 
and modulation efficiency of the AHWP depend on the relative angles between 
the plates in the stack and find that high modulation efficiency can be achieved 
with alignment accuracy of few degrees. 
© 2008 Optical Society of America 
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1. Introduction 

Recent experimental efforts in observational cosmology have been focused on searching for a 
signature from an inflationary period that occurred a short instant after the big bang. This 
signature is predicted to be imprinted in the polarization of the cosmic microwave background 
(CMB) radiation. Inflation predicts an inflationary gravitational- wave background (IGB) 
that left a particular pattern of linear polarization on the CMB. This pattern is different 
from that originating from primordial density anisotropy, which is the main source for the 
spatial intensity fluctuations over the sky and for a different, stronger linear polarization 
pattern. In the past few years, several groups have started to characterize the polarization 
signal coming from the primordial density anisotropy. 1 ' 2 The polarization pattern from the 
IGB is expected to be at least an order of magnitude smaller and it has not been detected 
yet. 

Thorough understanding of foregrounds and good control of systematic errors will be 
required to extract the small signal from the IGB. Both requirements lead to polarimeter 
designs that implement broad frequency coverage. Examples of such polarimeters are EBEX, 
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a NASA supported balloon-borne experiment 3 that is being constructed by the authors of 
this paper and their collaborators, BICEP/SPUD, CLOVER, QUIET, SPIDER, PAPPA, 
and PolarBear. 4-9 

A common technique to measure linearly polarized radiation is to use a rotating half-wave 
plate (HWP) together with a linear polarizer. The technique has been used extensively in 
the optical and IR wavelengths. 10-13 The first experiment to report CMB polarization results 
with this technique was MAXIPOL. 14,15 Although a HWP is a device that operates over 
narrow band of frequencies there are standard techniques to construct an 'achromatic HWP' 
(AHWP) that operates over a much broader range. 16-20 An AHWP is a stack of birefringent 
plates that are aligned with specific relative orientation angles between their optic axes. With 
an appropriate choice of angles it is possible to achieve modulation efficiency that is close to 
100 % over a large fraction of the millimeter wavelength band. 20 

The linear polarization content of incident radiation can be characterized in terms of two 
parameters, the degree of polarization P in and the orientation angle a in . An equivalent set 
is the normalized Stokes parameters Qm/hn and Ui n jli n (see Section [2] for a definition). To 
reconstruct these parameters from the signal detected by the instrument it is essential to 
quantify the extent to which the polarimeter itself changes the input parameters. 

In this paper, we present a computational study of the effects introduced by three- or 
five-stack AHWPs that are designed to fit CMB polarimeters operating in the range 120 — 
480 GHz. This work is motivated by EBEX that will implement a rotating 5-stack AHWP to 
search for the faint signals from the IGB. In Sections [2] and [3] we describe the mathematical 
formalism and define the figures of merit that are used to reconstruct the state of incident 
polarized light from the measured intensity. In Section [4] we use the figures of merit to 
quantify how well the incident polarization can be reconstructed. Section [5] discusses effects 
that arise from uncertainties in the spectrum of the incident radiation. In Section [6] we assess 
the performance of an AHWP as a function of its construction parameters. A summary of 
the key conclusions is given in Section [7j 

2. Polarimeter Model 

We consider a polarimeter that consists of an AHWP that rotates at a frequency /o, a 
linear polarizer, and a power detector (e.g. a bolometer), as shown in Figure [TJ The angle p 
gives the rotation angle of the AHWP around its axis of symmetry z. In such a polarimeter 
information about the incident polarization is contained in the intensity that is detected by 
the detector as a function of p. To a good approximation the detected intensity is sinusoidal 
as a function of p with a frequency of 4/ when there is a high signal-to-noise ratio. Our 
primary interest in this paper is to analyze the detected intensity as a function of p, which we 
call IVA (intensity vs. angle), with the purpose of reconstructing the incident polarization. 
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Fig. 1. A schematic diagram of the HWP polarimeter model. The transmission axis of a 
linear polarizer is parallel to the x axis. 



We use Mueller matrices to describe the signal in the approximation of normal incidence 
on the WP. We neglect the effect of absorption by the wave plate or effects of reflections 
between media that have different indices of refraction. 

Consider an input Stokes vector Si n of radiation propagating along the z axis that is 
incident on the polarimeter. The Stokes vector incident on the detector, which we call the 
output Stokes vector, is 

m 

Sam = Gl[[R{-p-6 i )r i {A6)R(p + 6 i ))S in {a in ,P in ), (I) 



i=l 

where 

A5 = 2ix-\n - n e \d, (2) 
c 

Sin Qini ^ ini 0) 

= , Pin COS sin2a:j n ,0), (3) 

G is the Mueller matrix of the linear polarizer, R is a rotation matrix, T is the Mueller matrix 
of a retarder, and Si n is the Stokes vector of the incident radiation, which is a function of the 
polarization angle aj n and of the degree of polarization P in . Information about the spectrum 
of the incident radiation is contained in I{v). We initially assume that the intensity of the 
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incident radiation is constant with frequency, I{y) = Iq = const . We discuss the effects of a 
non-constant incident spectrum in Section [5j Equation [T] assumes m wave plates in the stack; 
in this paper m — 1, 3 or 5. The variable AS is the retardance of a single wave plate and is a 
function of the ordinary and extraordinary indices of refraction n Q and n e , respectively, the 
thickness of a single wave plate d, and the electromagnetic frequency of light v. Throughout 
this paper, we assume that the incident radiation is not circularly polarized. We also assume 
that otin and P{ n are independent of v. The components of the Mueller matrices are 



° , (4) 
cos A5 — sin A5 
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As shown in Figure [TJ we choose the transmission axis of the ideal linear polarizer to be 
aligned with the +x axis. We define all the angles of rotation about the z axis with respect 
to the transmission axis of the grid. According to the usual convention, angles increase in 
the counter-clockwise direction from the +x axis in the xy plane. The relative orientation of 
plate i = 2, 3... in the stack relative to the first plate is given by 6\. The ordinary axis of the 
first plate is aligned with the x axis when p = degrees. We use the notation 6 to denote 
the entire set of relative orientation angles. 

The output of the detector is a function of its detection bandwidth and the intensity term of 
S ut- We assume that the detector has top-hat response of width Av about a center frequency 
v c . Note that the limit Av — > is equivalent to the case of illuminating the polarimeter with 
monochromatic light. With these assumptions the first element of the output Stokes vector 
can be written as 



(Iout){v c ,Av, 



2 



Av 
2 



I<mt{v, a in ,P in ,6,p)dv. 



(7) 
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A plot of {lout) as a function of p is the IVA. (Throughout this paper angle brackets ( ) 



and an analytic integration over any bandwidth is straight forward. For a three- and five- 
stack AHWP the analytic expressions are more complicated. 

It is the goal of this paper to discuss quantitatively how the amplitude and phase of 
the IVA depend on the construction parameters of the AHWP. Specifically, we make a 
quantitative mapping between the measured amplitude and phase of the IVA and the two 
parameters characterizing the incident polarization, the degree of input polarization P in and 
its orientation angle a.i n . 

When there is a finite detection bandwidth Av ^ the amplitude and phase of the IVA 
are calculated in the following way. We calculate the intensities as a function of angle p for 
each frequency within the bandwidth. We then sum the calculated intensities angle by angle 
to obtain a final IVA. The amplitude and phase are determined from that IVA. 

All of the analysis in this paper is computational. IVAs have been calculated as a func- 
tion of various parameters of the incident radiation and of the construction of the HWP. 
Many of the results were calculated by two independent computer codes to check for errors. 
Where practical, the results were compared to analytical calculations and agreement has 
been verified. 

3. Figures of Merit 

Figure [2] shows the IVA for the case of a single sapphire HWP and for an AHWP made of 
a stack of three and five sapphire plates. Table [T] lists the parameters of the plates used to 
generate these IVAs. The thickness of each wave plate gives AS = ir when vwp — 300 GHz. 
The top panels in Figure|2]show the IVA for a monochromatic detection bandwidth {Av = 0). 
In the panels on the bottom the detection bandwidth is Av = 60 GHz. The incident light is 
polarized parallel to the transmission axis of the grid, a in = degrees. 

Several generic features are apparent. The reduction in the amplitude of the IVA with 
a single HWP (left column) is a consequence of its chromaticity. Linear input polarization 
becomes elliptical when it passes through a wave plate that is optimized for a different 
frequency. There is substantially smaller reduction in amplitude of modulation for the 3-stack 
(middle column) or 5-stack (right column). However, whereas for a single plate the phase of 
the IVA is the same between different frequencies, or with a broad detection bandwidth, it 
becomes a function of frequency for the case of an AHWP. We define the phase angle <fi of 
the IVA as 
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{I out ) = A + A 4 cos (4p - 40), 



(9) 
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Fig. 2. IVA for monochromatic light (top panels) and for broadband radiation (bottom 
panels) for a single HWP, a three-stack AHWP, and a five-stack AHWP (left to right). See 
Table [l]for the parameters of the plates and for the details about the simulations used for the 
calculations. Frequencies of 150 (solid), 200 (dash), 250 (dot), 300 (dash-dot) GHz are used 
for the case of monochromatic light. For the broadband case we use 150 ± 30 GHz (solid) 
and 250 ± 30 GHz (dot). In all the panels, the maximum intensity is normalized to 1. 



where A and A 4 denote the average level and the modulation amplitude of the IVA, respec- 
tively. For the case of a single HWP the forms of A and A± are 



-TT o A5(u) 



A = ^ ^Av + P in cos 2a in J ^ cos 2 ' dv J , (10) 
and 

I Q P in p+¥ . 2 AS(v) 

Aa = / sin dv. 11 

2 2 v ' 



From a comparison of Equations [8] and [9] it is evident that in this case = <Xj n /2 and that 
it is independent of frequency. Since we chose a« n = degrees for the simulation shown 
in the left panels of Figure [2j <fi = degrees. However, for an AHWP the phase is a 
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Incident intensity 


/ = 1 


Indices of refraction of sapphire 21 


n Q = 3.047, n e = 3.364 


Thickness of each wave plate, d 


1.58 mm (<-> u WP = 300 GHz) 


Bandwidth of frequency, v c + Av 


150 ± 30 GHz, 250 ± 30 GHz 


Orientation angles of three-stack AHWP, #3 


(0, 58, 0) degrees 


Orientation angles of five-stack AHWP, #5 


(0, 29, 94.5, 29, 2) degrees 


Resolution of frequency 


0.5 GHz 


Resolution of wave plate angle 


0.1 degrees 



Table 1. Parameters of the wave plates and parameters used in the simulations to calculate 
the IVA. 



function of the thickness of the HWP, the detection bandwidth and the relative orientation 
angles. Mathematically <fi = 4>{v c , Au, cti n , 9), and therefore the IVAs in the middle and right 
columns of Figure [2] show non-zero phase angles. We define this overall 'phase offset' of the 
3- and 5-stack AHWPs as O - The quantity O is the value of when a in = degrees (e.g. 
0o ~ 65 degrees for the solid line of the middle bottom panel of Figure [2]). 

A useful figure of merit for the operation of a polarimeter is the 'modulation efficiency' 19 
defined as 

e = e(v c , Av,a in ,P in ,9) = (12) 

The efficiency e is a measure of the de-polarization introduced by the polarimeter and is an 
essential element in reconstructing the incident polarization Pj n from the one measured by 
the experiment P ou t- In this paper we calculate e by extracting P out from the IVA. P out is 
calculated from the ratio, 

J-, r-) 1 a r> n\ \-Lout) max {^out/min /-, \ 

Pout = Pontic, AU, OC in , P in , 9) = — r —— r . (13) 

\^out/max ~r \l-out/min 

Here (I ut)max and (I ut)min are the maximum and minimum of the IVA for angles < p < 90. 
(The modulation efficiency that was calculated in our earlier publication 20 assumed a some- 
what different functional form for P oni . See Section M for more details.) Using Equations [8j 



10 and 11 it is straight forward to show that for a single HWP and a single frequency 

Pi n sin , 

o A< 



2 A6(v) 

1 + P in cos 2a in cos^ — — 



and therefore 

p ■ 2 AS(u) 

/ \ lout Sill „ 

< U ) = -p— = V—p ~ ^A5M ■ ( 15 ) 

Pn 1 + P in cos 2a in cos z — ^ 



8 




Fig. 3. Modulation efficiency e = e(y, Al> = 0,cti n — 0,6) (top) and the phase offset 
0o = 4>{ain = 0, v, Av = 0, 8) (bottom) for the single HWP (left) and the three- (middle) 
and the five-stack (right) as a function of frequency. 



There are two cases for which this expression is particularly useful, (i) when P in is sufficiently 
small such that the denominator is approximately 1, and (ii) when «j n = 45 degrees (for any 
level of Pi n ). In both of these cases 

e{v) = sin (16) 

e is only a function of the retardance of the HWP and is independent of Pi n . In the first case 
e also does not depend on 

We note that instead of using Equation 13 a more generally appropriate process for ex- 
tracting the modulation efficiency is by fitting the IVA to a harmonic series of sine waves 
and then calculating 

e = (17) 

^0 
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Fig. 4. Top: Modulation efficiency of the single HWP, the three- and the five-stack AHWPs 
as a function of detection bandwidth for input polarization angle of (solid line), 22.5 (dot), 
45 (dash), 67.5 (dot-dash), and 90 (three-dot dash) degrees. Bottom: Output phase angle 
of the single, three-, and five-stack as a function of detection bandwidth for the same input 
polarization angles as the top panels. For both the modulation efficiency and the phase, 
v c = vwp- 



where Aq and A4 are the coefficients of the zeroth and fourth harmonic terms, respectively. 
For the results presented in this paper there is no difference between the two processes. 

The upper panels in Figure [3] give the modulation efficiency of a single HWP and of a 
three- and five-stack AHWP as a function of frequency. To calculate these efficiencies we 
analyzed monochromatic IVAs. The reduction in efficiency for the three-stack at a frequency 
of, for example, 150 GHz can be mapped to the smaller amplitude IVA in Figure [2] for the 
same frequency. The set of three panels shows that a larger number of plates in the stack 
gives a broader bandwidth of high modulation efficiency. 

In order to reconstruct the polarization angle oti n of the incident polarization it is essential 
to examine the phase angle of the IVA. The lower panels in Figure [3] show <fio as a function 
of frequency. They were also extracted from IVAs calculated with monochromatic detection 
bandwidths. The phase offset varies with frequency even near the center frequency v-wp- 

To gain some insight into the effects that we quantify in subsequent sections let us compare 
the middle panels of Figures 2 and 3 in more detail. The solid lines shown in Figure 2 
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were calculated for a center frequency of 150 GHz. The upper panel in Figure 2 shows the 
corresponding monochromatic IVA. The modulation efficiency of ~0.95 and phase angle of 
~70 degrees of that IVA can be read directly from the middle panels of Figures 3 at a 
frequency of 150 GHz. The lower panel of Figure 2 shows an IVA that has been calculated 
for a detection bandwidth of ± 30 GHz around 150 GHz. It was calculated by averaging the 
intensities at different frequencies. We used a frequency resolution of 0.5 GHz (see Table [I]). 
Each of these IVAs has a modulation efficiency and phase offset that can be read off from 
the middle panels of Figure 3. Both the modulation efficiency and the phase offset vary over 
the bandwidth. As a consequence, the resulting IVA is a superposition of sine waves with 
different amplitudes and phases. Hence the final IVA is also a sine wave, but its amplitude 
and phase depends on averaging intensities over frequencies. 

In the next Section we discuss how the efficiency e and the phase angle (f) depend on the 
construction parameters of the stack of HWPs, on the center frequency, and on the detection 
bandwidth, and how to relate them to the parameters of the incident polarization P in and 

4. Reconstruction of P in and aj n 

4- A. Modulation Efficiency and Phase 

The top panels of Figure [4] show the modulation efficiency of a single HWP, three-stack, and 
five-stack AHWPs as a function of bandwidth around v c = 300 GHz. The different curves 
correspond to different input polarization angles «j n . To calculate the modulation efficiency 
we used P« n = 1. For cti n = degrees modulation efficiency that is larger than 0.99 is achieved 
with a bandwidth of 200 (300) GHz for the three (five)-stack AHWP, while a single HWP 
achieves a bandwidth of only 50 GHz. For a given bandwidth the modulation efficiency is a 
function of the orientation of the incident polarization at n . So in order to reconstruct Pj n , 
information about ai n needs to be extracted first from the measured phase angle <f>. 

The bottom panels of Figure [4] show the output phase angle as a function of bandwidth 
around v c = 300 GHz. The phase of the single HWP shows flat response over the bandwidth. 
The phases of the IVA of the three- and five-stack AHWPs are a function of bandwidth, a 
result consistent with the bottom row of Figure [3j For a given bandwidth the phase angle 
has an overall offset 4>q. 

The conclusions so far are that if the incident radiation is known to be fully polarized and 
the detection bandwidth is known, then the orientation angle of the incident polarization 
and the modulation efficiency can be extracted. Alternatively, if the orientation angle of the 
incident fully polarized radiation is known then the modulation efficiency and an equivalent 
detection bandwidth can be extracted. These situations are encountered in the laboratory 
when calibrating the polarimeter. 
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Fig. 5. The extracted degree of polarization P out as a function of the degree of polarization 
of the incident light Pi n for the single-, three-, and five-stack. Each curve corresponds to 
the input polarization angle of (solid line), 22.5 (dot), 45 (dash), 67.5 (dot-dash), and 90 
(three-dot dash) degrees. The frequency and the bandwidth are v c ± Av = 150 ± 30 GHz 
(top) and 250 ± 30 GHz (bottom). For all the panels, v WP = 300 GHz. 



4-B. P out VS. Pi n 

The results in Figure [4] were calculated with the assumption of incident polarization Pj n = 1. 
We relax this assumption in Figure [5j which shows 

Pout 9* 

function of Pi n for various angles 
oti n . The local slope of each curve is the modulation efficiency e. The modulation efficiency 
is a function of both aj n and Pi n . That this is the case for a single HWP is evident from 
Equation 15 In an actual observation both Pj n and «j n are a-priori unknown, which suggests 
that reconstructing the polarimeter modulation efficiency, or the incident polarization Pj n is 
subject to additional uncertainty. In many practical cases this is not a limitation for reasons 
that we now discuss. 

Figure [5] shows that in cases where the incident polarization is known to be small the 
modulation efficiency is to a good approximation constant that does not depend either on 
Pi n or on a,in (for the simple case of a single HWP see Equations 15 and 16). The region 
where the approximation 'small incident polarization' is valid depends on the construction 
parameters of the HWP and the detection bandwidth. For example, for the top middle and 
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150 ± 30 GHz 


250 ± 30 GHz 


single HWP 




0.93l°;°ii 


three-stack AHWP 


0.89+gS 


0.996i°X 


five-stack AHWP 


0.9761°;°°? 


0.999™ 



Table 2. The modulation efficiency at P in — 0.1 with a.i n — 45 degrees is shown. The 
modulation efficiency is calculated as a slope of P out — P m relationship in Figure [5] The 
quoted errors are e mal - £45 and e m m — £45, where £45 corresponds to the modulation efficiency 
at ai n = 45 degrees at Pi n = 0.1. The maximum and the minimum modulation efficiency 
corresponds to ai n = 90 and degrees, respectively. 



bottom left panels the approximation is valid for P in < 0.2. It is valid for a much larger range 
of Pi n when using a 5-stack (see right panels), and even for the three-stack when it is used 
reasonably close to the designed band center (see middle bottom panel). 

We note that the HWPs that are used for the calculations shown in Figure [5] are each 
designed for a center frequency of 300 GHz (see Table [I]). Therefore the top left panel that 
shows the largest variation of the modulation efficiency with incidence angle is never likely 
to be used in practice. It describes a single HWP optimized for 300 GHz that is used for a 
band around 150 GHz. This panel is only shown for didactic purposes. 

We find then that in many practical situations there is a unique relation between P in and 
Pout, a relation that does not depend on the orientation angle of In the more general 
case when the value of P out depends both on P in and on ctj n , the value of ai n needs to be 
determined first from the IVA. This is straight forward for a single HWP because the phase 
<p of the IVA is equal to aj„/2 for any detection bandwidth; see for example the bottom left 
panel of Figure |4| The case of an AHWP is discussed in Section 4.C but the conclusion is 
that for a specified bandwidth there is a unique relationship between the phase and the 
angle a in . Therefore the procedure for finding P in is to first determine «j n using the phase 
of the IVA and then use the relation between P out and Pj n that is appropriate for this «j n . 

Laboratory measurements of modulation efficiency typically use incident polarizations 
that are close to Pj n = 1 in order to increase the signal to noise ratio of the measurement. 
Figure [5] demonstrates that determinations of e depend on the polarization angle «j n . An 
efficiency value that was determined in the laboratory using a particular angle «j n will not in 
general correspond to the modulation efficiency of the polarimeter during actual observations 
for which a^ n is not known. A simple remedy is to align the incident polarization in the 
laboratory such that «j n = 45 degrees. For that particular value the efficiency €45 is a 
constant as a function of Pi n and is equal to the same efficiency that would be measured 
with small incident polarizations. Table [2] summarizes this point in a quantitative way. The 
values shown give the efficiency expected with a in = 45 degrees for different frequency bands 
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Fig. 6. The output phase angle of the three- (left) and the five-stack (right) AHWPs as 
a function of the input polarization angle. The top panels give results for 150 (solid) and 
250 GHz (dot), each with a fixed bandwidth of ± 30 GHz. The bottom panels give results 
for a fixed center frequency of 300 GHz with bandwidths of ± (solid), ±100 (dot), and 
±200 (dash) GHz. 



and for different HWP configurations. The upper and lower values marked with ± give the 
additional increments of efficiency that would be determined if ai n was 90 (for plus) or (for 
minus) degrees. For example, the modulation efficiency of a single HWP (that is constructed 
according to the parameters given in Table [T]) at 150 GHz with a bandwidth of ±30 GHz is 
0.5 when measured with <Xj n = 45 degrees. This value of e does not depend on the magnitude 
of Pi n . Yet for observations with P in = 0.1 the modulation efficiency would be 0.055 (0.045) 
for «j n = 90 (0) degrees. Whereas for the single HWP the variation in modulation efficiency 
could be as large as 10%, it is about 2% or smaller with the three- or five-stack. These values 
depend on the construction parameters of the HWP and on the detection bandwidth and 
thus can not be taken as general. 

4.C. 4> vs. a in 



In Section 4.B we investigated how the measured degree of polarization P out relates to the 



input polarization P« n . We now quantify a similar relationship between and ai n . Figure [6] 
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shows the phase angle as a function of the input polarization angle for the three- and five- 
stack AHWPs. (Recall that <p is a constant over frequency for a single HWP, see the bottom 
left panel of Figure |1J) The panels show that a>i n and <p have a linear relationship with a 
slope of 0.5, and that this slope does not depend on the construction parameters of the HWP 
nor on the detection bandwidth. However, the phase offset is a function of v c and Av. 

The magnitude of the phase offset 0o is a critical parameter in the reconstruction of an 
unknown incident polarization angle Oj n . Since this phase offset is a function of the spectral 
response of the instrument it can either be calculated, if the spectral response is known, or 
measured in the laboratory by varying a iri of a known source and extracting from the IVA. 
(See however Section [5] for important caveats.) Errors in this calibration will propagate to 
errors in the determination of ctj n for a source whose polarization properties are not known. 

The direction of rotation of the AHWP affects the relationship between <p and oii n . With 
single HWP <fi = ±ctj n /2, where the sign is determined by the direction of rotation. For our 
particular choice of directions (see Figure [T]) we have <\> = +a in /2. However, the orientation 
angles of the stack of plates break the rotational symmetry for an AHWP and in general 
there are four possible choices 

= ±^«m + 0O, (18) 
1 7T 

= ±-a in + - - <p . (19) 



In Equation 19 the phase offset is no longer 0o> but | — 0q. For our particular choice, where 
both p and the AHWP orientation angles 9 are counterclockwise in the xy plane (as shown 
in Figure [T]), Equation 18 with a plus sign gives the relevant functional dependence. 



5. Spectrum of Incident Radiation 

So far we assumed an incident radiation spectrum that was constant with frequency. We 
now address the more general case where the spectrum of the incident radiation is a function 
of frequency. In this case the phase offset 0o depends on the details of this spectrum. Since 
0o is required for reconstruction of P, m and ctj n , the consequence is that knowledge of the 
incident spectrum is also required. 

To assess this effect quantitatively we consider three distinct spectra: (1) black body 
with the temperature of the cosmic microwave background radiation, (2) black body with 
a temperature of 300 K, and (3) galactic dust. We choose these spectra because they are 
relevant for calibration and for measurements of the polarization of the CMB at frequencies 
between 100 and 500 GHz. We assume the following spectra 

Icmb{v) = B(T CM B,v), (20) 
W") = A^B(T dust ,u), (21) 
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CMB 


Dust 


Lab 


150 ± 30 GHz 


0o = 57.86 


56.69 


57.33 


250 ± 30 GHz 


51.12 


51.16 


51.14 


420 ± 30 GHz 


53.85 


54.50 


54.49 





CMB - Dust 


CMB - Lab 


Dust — Lab 


150 ± 30 GHz 


A0 = 1.17 (Aa in = 2.34) 


0.53 (1.06) 


-0.64 (-1.28) 


250 ± 30 GHz 


-0.04 (-0.08) 


-0.02 (-0.04) 


0.02 (0.04) 


420 ± 30 GHz 


-0.65 (-1.3) 


-0.64 (-1.28) 


0.01 (0.02) 



Table 3. Top: The offset angles with four different spectra are shown. Bottom: The dif- 
ference of the offset phase between different spectra. The number in a parenthesis is the 
difference in terms of the polarization angle on n on the sky. A unit of the phase is in de- 
grees. 



I lah {v) = B(T lab ,u), (22) 

, m s 2irh z/ 3 . , 

B(T,v) = — -j- , (23) 

C e k B T _ I 

where B denotes a black body spectrum, T C mb = 2.73 K, A = 4 x 10~ 7 , 7 = 1.75, T dust = 
18 K, and = 300 K. We assume that the fractional polarization and the polarization 
angle of the incident radiation do not depend on frequency. We also assume that the degree 
of linear polarization of the CMB, of the galactic dust, and of a 300 K black body source are 
Pcmb = 1 x 10 -6 , Pdust = 0.1, and P iab = 1, respectively. 

The calculated phase offsets of a five-stack AHWP are summarized in Table [3j The top 
table shows the phase offset in unit of degrees. The bottom table shows the difference of the 
phase offsets between the different spectra. The parentheses indicate the level of difference 
in terms of polarization angle 

Assume that a 300 K source is used in the laboratory to calibrate the phase offset and that 
the laboratory measurement agrees with the phase offsets given in the right hand column 
of the top part of Table [3j If these values are used for either CMB or dust observations, 
they would give rise to errors in position angle of the polarization on the sky as given in 
parentheses in the two right columns of the bottom table. The correct prescription is to 
validate the design of the instrument using the laboratory measurements and then use the 
predicted phase offsets given assumptions or measurements of the spectra of the sources. An 
uncertainty in the knowledge of the spectrum would give an uncertainty in the determination 
of <Xi n . The designer of a polarimeter with an AHWP should plan for this uncertainty and 
its mitigation during the analysis of the data. 
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100 ZOO 300 400 500 100 ZOO 300 400 500 

Bandwidth [GHz] Bandwidth [GHz] 

Fig. 7. The modulation efficiency (left) and the phase offset (right) of the three-stack 
AHWP as a function of the angle of the second plate 02 and the bandwidth Av around a 
center frequency of vwp — 300 GHz. The color scale of the phase offset is in units of degrees. 
In both plots, the input polarization angle is ai n = 0. 




Bandwidth [GHz] Bandwidth [GHz] 

Fig. 8. The modulation efficiency (left) and the phase offset (right) of the five-stack AHWP 
as a function of the orientation angles of the second and fourth plates. The other angles are 
fixed at the values given in Table [TJ 
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Bandwidth [GHz] Bandwidth [GHz] 

Fig. 9. The modulation efficiency (left) and the phase offset (right) of the five-stack AHWP 
as a function of the orientation angle of the third plate. The other angles are fixed at the 
values given in Table [T] 




100 200 300 400 500 100 200 300 400 500 

Bandwidth [GHz] Bandwidth [GHz] 



Fig. 10. The modulation efficiency (left) and the phase offset (right) of the five-stack 
AHWP as a function of the orientation angle of the fifth plate. The other angles are fixed 
at the values given in Table [T] 
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6. AHWP Performance vs. orientation angles 9 

We have already pointed out in an earlier publication that it is relatively easy to achieve a 
high modulation efficiency with the three and five-stacks in terms of the requirements on the 
relative orientation of the plates. 20 In this section we expand on our earlier work and give a 
more thorough discussion. All of the analysis in this section assume a constant spectrum for 
the incident radiation. 

The left panels of Figures [7j [HJ |9j and [10] give contour plots for the modulation efficiency 
as a function of the orientation of the plates in the stacks. The modulation efficiency is 
calculated based on Equation 13 with Equations [7] and 12 The right panels of the same 
figures show the phase offset as defined in Equation |9j 

For a three-stack Title 18 showed that the highest modulation efficiency is achieved with 
a set of angle 9 = (0,58,0). The left panel of Figure [7] shows that this modulation is a 
weak function of the orientation of the middle plate near peak efficiency. The right panel of 
Figure [7] shows that for a second plate orientation 80 < # 2 ~ 100 degrees the phase offset O 
is essentially independent of detection bandwidth. This orientation angle, however, does not 
give the broadest range of frequencies for high modulation efficiency. On the other hand, with 
an angle of 58 degrees, which gives the broadest range of modulation efficiency, the phase 
offset has stronger dependence on the detection bandwidth. Experiment designers need to 
consider this trade-off between bandwidth for high modulation and for constant phase offset. 
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The three-stack has zero modulation efficiency at 6 2 close to 20 and 160 degrees and 
detection bandwidth of 400 GHz. This is because there is a strong variation of the phase 
offset angle <po with frequency near these parameters. Therefore, these points in the parameter 
space give the resultant IVA zero modulation amplitude, and correspondingly no phase can be 
defined as demonstrated by the singularities in the phase offset panel. The color discontinuity 
extending from the phase offset singularity toward bottom right is a consequence of phase 
offset periodicity. It is neither an artifact nor a real discontinuity. Phases that are larger 
than 90 degrees are interpreted as positive values close to zero. Similar features appear in 



Figures M M and 10 



With the five-stack, achieving high modulation efficiency requires higher accuracy of align- 
ment of the second, third, and fourth wave plates than that required in the case of the three- 
stack. Little accuracy is required from the orientation of the fifth plate in the five-stack. 
The efficiency is most sensitive to the orientation of the second and fourth plates, and an 
accuracy of 5 degrees is required to maintain efficiency higher than 0.95 over 300±150 GHz. 

In a previous publication 20 we gave results for the modulation efficiency that was based 
on the following expression 

P, _ ; I max Imin \ tOA\ 
out = \ J — j /, 1/4 J 

1-inax ~r -'mm 



which is different from the more correct definition given in Equation 13 The left panel 



of Figure [TT] is the modulation efficiency based on Equation [24] with the same parameters 
that produced Figure [7| The right panel shows the differences between the two results. The 
modulation efficiency in Figure [7] accounts for the phase variation of the IVA curves as a 



function of frequency. In contrast the modulation efficiency in Figure 11 does not encode 
this variation. 



7. Discussion 

We analyzed the performance of a three- and five-stack AHWP polarimeters operating in 
the sub-millimeter wave band. Let us summarize the points that have been discussed and 
make some additional comments where appropriate. 

• Three- and five-stack AHWP polarimeters provide broad bandwidth with high modu- 
lation efficiency. 

• Their IVA has a phase offset that depends on the construction parameters of the stack, 
on the spectral response of the instrument, and on the spectrum of incident radiation, 
(our discussion assumed that the degree of polarization and the angle were independent 
of frequency with the detection bandwidth.) 

• If the spectral response of the instrument, and the spectrum of the source are known, 
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then measurements of the phase of the IVA can give the orientation angle of the incident 
polarization. 

We note that in many cases much of the radiation incident on the detector is due to 
emission by the telescope itself. If this emission is polarized it too will affect the phase 
of the IVA and hence the measurement of the angle of incident polarization. 

• measurement uncertainties in either the spectral response of the instrument or the 
spectrum of the source translate to uncertainties in the reconstruction of the angle of 
incident polarization. The amount of uncertainty needs to be assessed on a case-by-case 
basis. 

• Measurements of the amplitude of the IVA, which gives the degree of output polar- 
ization, can be uniquely inverted, in most cases, to give the input polarization if the 
modulation efficiency is known. 

• In some cases information about the angle of the incident polarization needs to be used 
together with the modulation efficiency to find the incident degree of polarization. 

• Laboratory measurement to find the modulation efficiency that are conducted with a 
source that has high degree of polarization should have an incident polarization angle 
of 45 degrees. At this angle the measured efficiency is the same as would be measured 
at any angle when Pj n is small. 

• We discussed how the modulation efficiency and phase offset of the polarimeters depend 
on errors in the orientation of the plates. Generally, an accuracy of few degrees is 
sufficient to ensure close the ideal performance. 

• We discussed the how the incident spectrum of the radiation affects the IVA and the 
extraction of the parameters of the incident radiation. 

Our analysis assumed a spectral response of the instrument that was top-hat in shape 
over a range in frequencies. This is an idealization. In any practical instrument, the 
entire spectral response of the instrument is necessary in order to reconstruct the 
parameters of the incident polarization. 
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List of Figure Captions 

Fig. 1. A schematic diagram of the HWP polarimeter model. The transmission axis of a 
linear polarizer is parallel to the x axis. 

Fig. 2. IVA for monochromatic light (top panels) and for broadband radiation (bottom 
panels) for a single HWP, a three-stack AHWP, and a five-stack AHWP (left to right). See 
Table [T] for the parameters of the plates and for the details about the simulations used for the 
calculations. Frequencies of 150 (solid), 200 (dash), 250 (dot), 300 (dash-dot) GHz are used 
for the case of monochromatic light. For the broadband case we use 150 ± 30 GHz (solid) 
and 250 ± 30 GHz (dot). In all the panels, the maximum intensity is normalized to 1. 

Fig. 3. Modulation efficiency e = e(u, Av = 0, a in = 0, 0) (top) and the phase offset 
0o = 0(«m = 0, is, Av = 0, 6) (bottom) for the single HWP (left) and the three- (middle) 
and the five-stack (right) function of frequency. 

Fig. 4. Top: Modulation efficiency of the single HWP, the three- and the five-stack AHWPs 
as a function of detection bandwidth for input polarization angle of (solid line), 22.5 (dot), 
45 (dash), 67.5 (dot-dash), and 90 (three-dot dash) degrees. Bottom: Output phase angle 
of the single, three-, and five-stack as a function of detection bandwidth for the same input 
polarization angles as the top panels. For both the modulation efficiency and the phase, 
v c = vwp- 

Fig. 5. The extracted degree of polarization P out as a function of the degree of polarization 
of the incident light Pj ra for the single-, three-, and five-stack. Each curve corresponds to 
the input polarization angle of (solid line), 22.5 (dot), 45 (dash), 67.5 (dot-dash), and 90 
(three-dot dash) degrees. The frequency and the bandwidth are v c ± Av = 150 ± 30 GHz 
(top) and 250 ± 30 GHz (bottom). For all the panels, v-wp — 300 GHz. 

Fig. 6. The output phase angle of the three- (left) and the five-stack (right) AHWPs as 
a function of the input polarization angle. The top panels give results for 150 (solid) and 
250 GHz (dot), each with a fixed bandwidth of ± 30 GHz. The bottom panels give results 
for a fixed center frequency of 300 GHz with bandwidths of ± (solid), ±100 (dot), and 
±200 (dash) GHz. 

Fig. 7. The modulation efficiency (left) and the phase offset (right) of the three-stack 
AHWP as a function of the angle of the second plate 82 and the bandwidth Au around a 
center frequency of v WP = 300 GHz. The color scale of the phase offset is in units of degrees. 
In both plots, the input polarization angle is a in = 0. 

Fig. 8. The modulation efficiency (left) and the phase offset (right) of the five-stack AHWP 
as a function of the orientation angles of the second and fourth plates. The other angles are 
fixed at the values given in Table [Tj 

Fig. 9. The modulation efficiency (left) and the phase offset (right) of the five-stack AHWP 
as a function of the orientation angle of the third plate. The other angles are fixed at the 



24 



values given in Table [T} 

Fig. 10. The modulation efficiency (left) and the phase offset (right) of the five-stack 
AHWP as a function of the orientation angle of the fifth plate. The other angles are fixed at 
the values given in Table [1} 



Fig. 11. The modulation efficiency of a three-stack AHWP based on Equation 24 (left) 



and the difference between this efficiency and the one calculated in Figure. 11 (right) 



Table Caption 

Tab. 1. Parameters of the wave plates and parameters used in the simulations to calculate 
the IVA. 

Tab. 2. The modulation efficiency at P« n = 0.1 with oti n = 45 degrees is shown. The 
modulation efficiency is calculated as a slope of P out — Pi n relationship in Figure [5] The 
quoted errors are e max — e 45 and e min — e 45 , where e 45 corresponds to the modulation efficiency 
at ai n = 45 degrees at Pj n = 0.1. The maximum and the minimum modulation efficiency 
corresponds to = 90 and degrees, respectively. 

Tab. 3. Top: The offset angles with four different spectra are shown. Bottom: The difference 
of the offset phase between different spectra. The number in a parenthesis is the difference 
in terms of the polarization angle a« n on the sky. A unit of the phase is in degrees. 
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